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The ground state of a bilayer quantum Hall system at i*r = 1 with model pseudopotential 
is investigated by the DMRG method. Firstly, pseudopotential parameters appropriate for the 
system with finite layer thickness are derived, and it is found that the finite thickness makes the 
excitonic phase more stable. Secondly, a model, where only a few pseudopotentials with small 
relative angular momentum have finite values, is studied, and it is clarified how the excitonic 
phase is destroyed as intra-layer pseudopotential becomes larger. The importance of the intra- 
layer repulsive interaction at distance twice of the magnetic length for the destruction of the 
excitonic phase is found. 

KEYWORDS: bilayer quantum Hall system, excitonic state, phase transition, composite fermion, strong 
magnetic field, DMRG 



1. Introduction 

The bilayer quantum Hall system at total filling fac- 
tor zAp = 1 is known to be in the inter-layer coherent 
state when the layer separation is small enough. 1-5 ) The 
coherent state in the limit of vanishing layer separation 
is described by Halperin's state, 6 ' which can also 

be considered as an excitonic state, where an electron in 
one of the layer is bound to a hole in the other layer. 
Many experimental evidences for this excitonic correla- 
tion have been accumulated. 7_11 ) The excitonic correla- 
tion becomes gradually smaller as the layer separation 
increases. It has been clarified that the excitonic correla- 
tion almost vanishes at around (d/l) c = 1.85 for typical 
samples, 7 ~ 10 ) where d is the distance between the two lay- 
ers and I = \Jh/eB is the magnetic length. 12 ) At larger 
layer separation, the two layers become independent of 
the other layer, and each layer is in the composite fermion 
(CF) liquid state. 13 - 15 ) 

The details of the transition between the excitonic 
phase and the independent CF liquid state, however, 
have not been clarified in spite of intense theoretical 16-27 ) 
and experimental investigations. 28 ~ 34 ) There are several 
issues to be clarified both theoretically and experimen- 
tally. The first issue is what is the critical value of d/l 
at which the excitonic phase is terminated. The second 
issue is whether the transition is continuous or discontin- 
uous. The third issue is whether electron spin is involved 
in this transition or not. Actually, these issues are related 
to each other. 

As for the involvement of the spin degree of freedom, 
all the theoretical investigations have assumed full spin 
polarization both in the CF liquid phase and in the ex- 
citonic phase. This is because the system becomes much 
simpler if we can assume full spin polarization, and we 
can expect the spin degree of freedom can be suppressed 



in the strong magnetic field where Zeeman splitting is 
large. Furthermore, full spin polarization in the exci- 
tonic phase is expected even in smaller Zeeman splitting 
from the following theoretical consideration. Namely, it 
should be noticed that the excitonic state can also be 
understood as pseudospin ferromagnetic state, where the 
z component of the pseudospin stands for the layer in- 
dex. 1 ) The Coulomb interaction between electrons in the 
different layers forces the pseudospin align in the xy- 
plane, and the excitonic state is realized. Considering 
the facts that the interaction between electrons with op- 
posite spin is stronger than that between electrons in 
different layers at finite d/l, and the existence of the 
Zeeman splitting, real spin of the electron must also be 
ferromagnetically aligned when the excitonic state (i.e. 
pseudospin ferromagnetic state) is realized. Thus, the 
theoretical estimates for the boundary value (d/l) c have 
been done between spin-polarized excitonic phase and 
the spin-polarized CF liquid phase. This assumption also 
leaves a possibility for continuous transition between the 
two phases. 

With this assumption of the full spin polarization, the 
estimates of (d/l) c is not so different from theory to the- 
ory. However, even in this case, there has been no con- 
sensus on the details of the transition. Some of the the- 
ories suggest a first order phase transition, 20 ' 22 ' 23 ) and 
the others suggest second order phase transition. 26 ) The 
possibility of smooth crossover between the two limiting 
states has not been denied. 25,27 ) Possibility of interme- 
diate states is also raised. 16 ~ 18 ) 

The situation become quite different if the CF liquid 
phase is partially spin-polarized. The transition between 
the two states with different spin polarization cannot be 
continuous. Moreover, the estimate of d/l should be re- 
considered. This is because the partial spin polarization 
in the CF liquid phase means that the spin-polarized CF 
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liquid has higher energy, so the boundary between the 
excitonic phase should have moved to the lower value of 
d/l than the estimate assuming full polarization. 

Now, recent experiments are giving evidences that the 
spin degree of freedom is important. 29 ' 30 ' 32,33 ) It has 
been found that there is a first order phase transition 
between excitonic state and partially spin-polarized CF 
liquid state in typical experimental situations, although 
the transition looks like continuous possibly due to in- 
homogeneity of the sample. This fact has raised three 
questions: (i) what property the state in the partial spin 
polarization has? (ii) how the transition occurs if the Zee- 
man splitting is enhanced? and (iii) what is the value of 
(d/l) c , if enhanced Zeeman splitting makes both phase 
spin-polarized? For the first question theoretical investi- 
gation has not been developed, mainly because the na- 
ture of the partially spin-polarized CF liquid state has 
not been clarified. For the second and third questions, it 
has been found experimentally that the excitonic phase 
persists to larger layer separation, when the spin Zeeman 
splitting is enhanced by in-plane magnetic field, and the 
transition between polarized states now becomes contin- 
uous. 32,33 ' In one of the experiments, value as large as 
(d/l) c ~ 2.3 is reported. 32 ' 

The large value of the transition point for spin- 
polarized system is surprising, because theoretical es- 
timate for the transition point is around d/l = 1.8 
or smaller. 5,16,20,26,27 ^ We previously have investigated 
fully spin-polarized system by DMRG method. 26 ) In that 
paper we considered a system where the thickness of each 
layer is negligibly small. We calculated the excitonic cor- 
relation, and found that it decreases as the layer sep- 
aration becomes larger, and almost vanishes at around 
d/l = 1.6. 

In the present paper, we restrict ourselves to fully spin- 
polarized systems as before. We try to understand the 
shift of the phase boundary to larger value of d/l, and 
try to understand phase transition in such spin-polarized 
systems. We suspect this discrepancy of the transition 
point between theory and experiment comes from the 
thickness of the two-dimensional layer. Thus, in section 
2 we calculate pseudopotential parameters taking into 
account thickness effect, and using those parameters we 
determine the transition point where the excitonic cor- 
relation vanishes. In order to investigate the mechanism 
of the transition in the spin-polarized system further, we 
consider in section 3 a minimal model in which only a 
few pseudopotential parameters are included. Discussion 
and summary of the results are given in section 4. Inves- 
tigation of the partially spin-polarized CF liquid state is 
left for the future work. 

2. Effect of Layer Thickness 

2.1 Pseudo-potentials 

Previously, development of the ground state as the 
layer separation increases has been mostly studied ne- 
glecting the thickness of the quantum wells making up 
the bilayer system. 35 ) In this paper we calculate Hal- 
dane's pseudopotentials 36 ) considering the layer thick- 



ness, and investigate the effect on the ground state. 

We consider a model bilayer system where each layer is 
realized as a quantum well with infinitely high confining 
potential. The width of the well is denoted as w, and the 
separation between the layers d is measured between the 
center of each well. We take the z-axis perpendicular to 
the 2-d layers, so that the electrons can reside either at 
-(d + w)/2 < z < -{d - w)/2, or at (d - w)/2 < z < 
(d+w)/2 in this model. The single particle wave function 
is written as 



®±{x,y,z) = cp(x,y)ip±(z), 
where + and — distinguish the layers and 
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Using these wave functions, we calculate Fourier trans- 
form of the inter-layer and intra-layer Coulomb interac- 
tion potentials, V e (q) and V a {q). 
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where e is the dielectric constant, a and a' are either + 
or — , a = a' for intra-layer interaction V a (q), and a ^ a 1 
for inter- layer interaction V e (q) . 

The integral can be done analytically. The inter-layer 
interaction V c (q) is calculated as follows 
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On the other hand, the intra-layer interaction is cal- 
culated as follows 37,38 ) 
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Haldane's pseudopotential is obtained from 

V*( g ) : 3«) 
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This integral is evaluated numerically. The results are 
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Fig. 1. Relative deviation of the inter-layer pseudopotcntials as 
functions of w/l for a system with d = 21. AV^ = (V r \\/V r7 \(w = 
0) — 1) X 100 are plotted as a solid line, dashed line, dash-dotted 
line, dottcd-linc, and dash-dot-dotted line for m = 0, 1, 2, 3, 4, 
respectively. 




w/l 

Fig. 2. Relative deviation of the intra-layer pseudopotcntials as 
functions of w/l. AV£ = (V£/V£(w = 0) - 1) x 100 are plotted 
as a solid line, dashed line, and dash-dotted line for m = 1,3, 5, 
respectively. 

plotted in Figs. 1 and 2 as relative deviation from the 
value at w = in percent as functions of w/l. Namely, 
AV£ = {V 7 c n /V 7 c n {w = 0) - 1) x 100 for m = 0, 1, 2, 3, 4 are 
plotted in Fig. 1, and AV^ for m — 1, 3, 5 are plotted in 
Fig. 2. We can see that the effect of finite w is larger for 
smaller relative angular momentum to. The inter-layer 
pseudopotentials increases for m = 0, 1, 2, but other 
pseudopotentials decrease as w increases. 

2.2 Results for the excitonic correlation 

We obtain the ground state of the bilayer system of 
finite electron number by DMRG method. 39 ~ 41 ) In this 
section we diagonalize rectangular systems with periodic 
boundary conditions, namely systems on torus, at to- 



Fig. 3. Excitonic correlation for double quantum well of layer 
thickness w = and w/d = 0.64. 

tal filling factor = 1. In order to see how the layer 
thickness affects the stability of the excitonic phase, we 
calculate excitonic correlation function for a system with 
finite w, and compare the results with our previous re- 
sult obtained for a system with w = 0. 26 > The excitonic 
correlation is defined by 

9UX) = N % N ^ ^\ c U^xcl afl c^), (7) 

where is the ground state and x is the creation 
operator of the electron in the a layer at center coor- 
dinate X, N is the total number of electrons, which 
is also the total number of single electron states per 
layer, N\ and N2 are number of electrons in each layer 
(Ni =N 2 = N/2). 

In Fig. 3, we show results of the excitonic correlation 
g ex (L x /2) for systems with w = and w/d = 0.64, where 
L x is the liner dimension of the system with periodic 
boundary condition. The total number of electrons N is 
20. The value of w/d is chosen to reproduce typical ex- 
perimental situation where d — 28 nm and w = 18 nm. 
As seen from the figure, the value of d/l where the ex- 
citonic correlation almost vanishes changes from 1.8 for 
w = to 2.0 for w/d = 0.64. 

The reason for the increase of the excitonic correla- 
tion and shift of the phase boundary as w increases is 
understood from the behavior of the pseudopotentials. 
The increase of inter-layer pseudopotcntial V^ is favor- 
able for the realization of the excitonic phase, because 
of the fact that the excitonic phase is the exact ground 
state for system with Vq > 0, and all other pseudopo- 
tentials being zero. The decrease of V* is also favorable 
for the excitonic phase, since this interaction is respon- 
sible for the destruction of the excitonic phase as we will 
see in the next section. The effect of finite w through 
other pseudopotentials should be small, since the devia- 
tion from the values at w = is much smaller for these 
pseudopotentials compared with Vq and V* 
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2.3 Comparison with experiment 

First, we examine how samples with different values of 
w/l give different critical values (d/l) c . We compare ex- 
periments where the Zeeman splitting is not enhanced 
artificially. Thus the transition may be between spin- 
polarized excitonic phase and partially spin-polarized CF 
liquid state, but qualitative effect can be seen. The sam- 
ple with smallest w/d is used by Wiersma et al. 28 ^ For 
this sample with w/d = 0.58, the critical value (d/l) c is 
1.65, and here w/l = 0.96. On the other hand, for typical 
samples with w/d = 0.64, the critical value (d/l) c is 1.85, 
and w/l = 1.18. These experimental results qualitatively 
agree with the present theoretical result. 

For more quantitative comparison, we need to compare 
theoretical result with experiments where the Zeeman 
splitting is enhanced by in-plane magnetic field such that 
both phases are completely spin-polarized. There have 
been two such experiments. Giudici ct al. 32 ) measured 
the activation energy of the quantum Hall state, which is 
identified as the excitonic phase, at various values of d/l. 
They found that as the Zeeman splitting is enhanced, the 
boundary for the quantum Hall state shifts from (d/l) c ~ 
1.9 at no enhancement to 2.3 for sufficiently enhanced 
Zeeman splitting. This value of the phase boundary is 
considerably larger than our result of (d/l) c = 2.0. 

On the other hand, in a more recent experiment, Finck 
et al. 33 ) investigated the phase boundary using Coulomb 
drag experiment. They also find increase of (d/l) c as 
in-plane magnetic field is increased, but the value of 
{d/l) c by their definition starts from 1.75 and saturates 
at around 1.85. This value is a little lower than our result. 
However, Fig. 1 of rcf. 33 allows us another way to de- 
termine the critical value. Namely, we can use the value 
of d/l at which the Hall drag resistance begins to have fi- 
nite value. Considering the fact that the Hall drag should 
be caused by finite excitonic correlation, this definition 
will be closer with our definition of (d/l) c . Adopting this 
definition we find that the critical value increases from 
(d/l) c — 1-8 in the perpendicular magnetic field to 2.0 in 
the magnetic field tilted by 66°. In this case agreement 
with our result is quite good. 

The samples these authors used have almost identical 
values of d and w. We do not know the origin of the dif- 
ference. It may be due to slight difference in the sample 
quality and/or geometry, or may be due to the differ- 
ent methods and definitions to determine the transition 
point. It has been pointed out that the values of d and 
w may not be exact. 29 ' 32 ) Anyway, we suspect that the 
most important factor making the difference in (d/l) c 
between the theory and experiment if any is an addi- 
tional effect caused by tilting of the magnetic field. In the 
experiments, the application of the in-plane field makes 
the total magnetic field tilting about 60 degrees from 
the direction perpendicular to the 2-d plane. The single- 
electron wave function tends to align with the magnetic 
field line. Namely, the wave function no longer has the 
separable form of eq.(l), and looks as if the center coor- 
dinate X shifts in the direction of the in-plane field as z 
increases within each quantum well. 
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Fig. 4. Excitonic correlation for Vq-Vi model in torus geometry. 
The total number of electrons TV is 18, and the aspect ratio 
Lx/Ly is 1. 

This deformation of the wave function brings 
anisotropy in the 2-d plane, and relative angular momen- 
tum no longer is a good quantum number. Investigation 
of such a case is not considered in this paper, although 
part of the effect may be able to be taken into account 
by calculating angle averaged effective pseudopotentials. 

3. V Q e -V* Model 

In this section we investigate a minimal model that 
reproduces phase transition between the excitonic phase 
and the CF liquid phase. The excitonic phase is realized 
as the exact ground state for a model where only the 
inter-layer short-range interaction exists between elec- 
trons. Namely, it is a model where only V e has finite 
positive value, and other pseudopotentials are zero. On 
the other hand, the short-range Coulomb repulsion is es- 
sential for the realization of the CF liquid state in a single 
layer at half-filling. Namely, a model where only has 
finite positive value, and other pseudopotentials are zero 
has ground state where two layers are independent, and 
in the CF liquid state. Therefore, a model in which only 
Vq and V* have positive values is the minimal model 
that has both phases as limiting cases. We characterize 
this model by a single parameter Vfi/V*. 

The systems we consider in this section are a finite 
size system on a torus (a rectangular system with peri- 
odic boundary conditions) and a system on a sphere. The 
results of both systems should be the same in the limit 
of large system. However, in finite systems, size effect 
appears through the boundary conditions. In rectangu- 
lar systems, the boundary conditions are specified by the 
area of the unit cell L x L y and the aspect ratio L x /L y , 
whereas spherical geometry is characterized only by the 
radius of the sphere R. We therefore remove the size ef- 
fect by a simple extrapolation of R in spherical geometry. 

The results for the excitonic correlation and charge gap 
are shown in Figs. 4-7. In Fig. 4, excitonic correlation for 
a square system with 18 electrons is plotted as a function 
of Vq/Vi- In the limit of Vfi /V* — > oo, the excitonic cor- 
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Fig. 5. Excitonic correlation for Vq-V^ model in torus geometry. 
N is the total number of electrons. 

relation <7 ex (L x /2) is unity. The excitonic correlation de- 
creases with the decrease in Vq/V*, and almost vanishes 
at around V e /Vf ~ 0.4. Similar behavior is obtained for 
other systems with different N and L x /L y . The decrease 
in the excitonic correlation is continuous and the ground 
state changes to the CF liquid state at some finite value 
of V£/Vf. 

Figure 5 shows a close-up view of the excitonic corre- 
lation near the transition point. In this figure, those for 
several other systems are also plotted to see size depen- 
dence. In the region of large Vq/V*, g ex (L x /2) decreases 
with the increase in L x /L y . Although excitonic correla- 
tion is independent of the distance between the two elec- 
trons in the limit of Vq/Vi = °°j feW has r dependence 
for finite value of V e /V® and g ex (L x /2) decreases with 
the increase in the distance L x /2. Around the transition 
point, however, different size dependence appears. In this 
region Vq/V* is small and the intra-layer interaction, V 1 a , 
is relatively important. Since the interaction energy from 
Vi is sensitive to the spatial configuration of electrons in 
the unit cell, the ground state energy of finite system 
depends on boundary conditions L x /L y . For the system 
of 16 electrons, each layer has 8 electrons. The energy 
from Vi has a minimum at L x /L y ~ 2 (2 x 4 configura- 
tion). For 18 electrons, each layer has 9 electrons. Then 
a minimum of the energy from appears at L x /L y ~ 1 
(3x3 configuration). Such /Independence enhances 
or reduces the effect of V* and the transition (crossover) 
point from the excitonic phase to the CF liquid phase is 
modified depending on N and L x /L y . 

To remove such size dependence, we calculate g e x(r) 
in spherical systems and extrapolate the results to the 
limit of large distance r and radius R as shown in the 
inset of Fig. 6. The extrapolated results are plotted in 
the main panel of Fig. 6. Here, g e x(r = oo) is the ex- 
citonic correlation between the electrons whose distance 
r is infinity, that means g ex {r = oo) represents macro- 
scopic coherence. Similarly to g cx (L x /2) of finite sys- 
tems, macroscopic coherence continuously decreases with 
the decrease in Vq/V^ and almost vanishes at around 



Fig. 6. Excitonic correlation for Vq-V* model in spherical geome- 
try. <?ex(»" = oo) is the extrapolated value of excitonic correlation 
in the limit of r — ¥ oo. Inset shows <? e xM for several different N. 
Dots in the inset for N = 24 and Vq /V* = 5 represent results 
in torus geometry of L x /L y = 1. Other results in the inset are 
calculated in spherical geometry. 

V e /V? ~ 0.4. 

Finally in Fig. 7 charge excitation gap A c is plotted. 
In the limit of Vq = oo, infinite energy is required to 
add or remove an electron, because at least one exciton 
consisting electron-hole pair is destroyed. This charge ex- 
citation energy is proportional to V e , and A c decreases 
with the decrease in V e . At the transition to the CF liq- 
uid state, A c should vanish in the limit of large system. 
In finite systems, however, discrete energy levels leave fi- 
nite value of charge gap even in the CF liquid state. In 
such a case, the charge gap has considerable size depen- 
dence and it tends to decrease with the increase in the 
system size. The charge gap in Fig. 7 shows large size 
dependence for small Vq/V^ indicating that the system 
is in the CF liquid state, but such size dependence is re- 
moved for larger Vq/V*. The change in the behavior of 
the size effect occurs at around V e /Vi ~ 0.4 and the gap 
seems to close in large systems for V e /Vi < 0.4. We also 
notice that between 0.4 < Vq/V* < 1.0 the charge gap 
deviates from the linear behavior at 1.0 < Vq/V^. This 
indicates that the excitonic correlation becomes weaker 
in this regime, and the energy of a charge excitation such 
as meron pair 5 -* that does not destroy an exciton be- 
comes lower than that of a simple charge excitation that 
destroys an exciton. 

These results clearly show that the excitonic correla- 
tion characterizing the excitonic phase decreases with the 
decrease in Vq/V^, and continuously vanishes at around 
Vq/Vi — 0-4 with the transition to the CF liquid state. 
This behavior is quite similar to the phase transition 
(crossover) of the actual system with Coulomb interac- 
tion as the separation becomes larger. 

4. Discussion 

We have seen that the excitonic phase stabilized by 
the inter-layer pseudopotential V C is destroyed by intra- 
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Fig. 7. Charge gap for Vq-V* model in spherical geometry. N is 
the total number of electrons. 



layer pseudopotential V*. In the excitonic phase, every 
single electron state is occupied that is constructed as a 
superposition of the wave functions of the same center 
coordinate from two layers: 

1 



-n>: 



x 



t 



+ e l °c 



-a,X 



)|0), 



(8) 



where is an arbitrary phase, and X is the center co- 
ordinate in the Landau gauge. The same wave function 
can be written in symmetric gauge, in which the angular 
momentum is a good quantum number. In this ground 
state, for an electron in a single electron state of cr-layer 
at quantum number X, there is no electron in the other 
layer at the same X. Namely, an electron at X in cr-layer 
is always bound to a hole in the —cr-layer. Thus, this 
is a state where the pseudopotential Vq is completely 
avoided. On the other hand, an electron has a nearest 
neighbor electron in the same layer with relative angular 
momentum 1 with probability 1/2 in this state. Thus, 
the pseudopotential makes the ground state energy 
higher. In the CF liquid state, this interaction is partly 
avoided as each electron bounds a hole to form an in- 
plane dipole. Namely, if we compare the pair correlation 
function between the excitonic phase and the CF liquid 
phase, the CF liquid phase has smaller value at around 
r/l = 2 where the pseudopotential has a peak. 26 ) The 
intra-layer pseudopotential V* acts to move the hole in 
the opposite layer to the same layer, and weakens the 
excitonic correlation. This mechanism of the phase tran- 
sition is what we suggested in the previous paper 
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the essence of the transition is reproduced in the present 
minimal model. 

We have investigated models where Vf and are 
included also. We have found that the excitonic phase 
is stabilized by adding positive Vf, but destabilized by 
adding positive V 5 a . This behavior is understood consid- 
ering the mechanism of the phase transition and the 
shape of the pseudopotentials in the real space. The 
two-electron state with relative angular momentum one 
(m = 1) has a peak probability for a configuration in 



which two electrons are separated by distance r = 21. 
The potential V* has a peak at this separation, but Vf 
has a dip here. 42 ) Thus, has an effect to weaken the 
repulsion of at r = 21, and reduces the pressure to 
bound a hole at this distance. On the other hand, V$ has 
a peak around r = 21, so positive V 5 a is harmful to the 
excitonic phase. These results show importance of the 
intra-layer interaction at distance r = 21 for destruction 
of the excitonic phase. The excitonic phase cannot have 
a hole in the same layer at distance r = 21, so the repul- 
sive interaction at this distance destroys the excitonic 
correlation. 

In §2, we have seen that the layer thickness stabilizes 
the excitonic phase to larger value oid/l. This is a reason- 
able result, since main consequence of the layer thickness 
is to make V e larger and to make V* smaller. 

In summary we have shown that the layer thickness 
has an effect to shift the phase boundary to larger value 
oid/l. It qualitatively explains the experimental fact that 
sample with small (large) w/l has smaller (larger) (d/l) c . 
For typical sample with w/l = 0.64, quantitative compar- 
ison is made. We partly explained (d/l) c larger than 1.8 
observed in recent experiments 32, 33 - ) in which the Zee- 
man splitting is enhanced by in-plane magnetic field. We 
also investigated a minimal model that show the transi- 
tion between the excitonic phase and the CF liquid state, 
and clarified importance of V* or the repulsive interac- 
tion at distance r ~ 21 for destruction of the excitonic 
phase. 
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